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ABSTRACT 

We present a study of the structural and scaling properties of the gas distributions in the intracluster medium (ICM) of 31 nearby 
(z < 0.2) clusters observed with XMM-Newton, which together comprise the Representative XMM-Newton Cluster Structure Survey 
(REXCESS). In contrast to previous studies, this sample is unbiased with respect to X-ray surface brightness and cluster dynamical 
state, and it fully samples the cluster X-ray luminosity function. The clusters cover a temperature range of 2.0 - 8.5 keV and possess 
a variety of morphologies. The sampling strategy allows us to compare clusters with a wide range of central cooling times on an 
equal footing. We applied a recently developed technique for the deprojection and PSF-deconvolution of X-ray surface brightness 
profiles to obtain non-parametric gas-density profiles out to distances ranging between 0.8i?joo and 1. 5^^500- We scaled the gas density 
distributions to allow for the systems' differing masses and redshifts. The central gas densities differ greatly from system to system, 
with no clear correlation with system temperature. At intermediate radii (~ O.SRsoo), the scaled density profiles show much less 
scatter, with a clear dependence on system temperature. We find that the density at this radius scales proportionally to the square 
root of temperature, consistent with the presence of an entropy excess as suggested in previous literature. However, at larger scaled 
radii this dependence becomes weaker: clusters with kT > 3 keV scale self-similarly, with no temperature dependence of gas-density 
normalisation. The REXCESS sample allows us to investigate the correlations between cluster properties and dynamical state. We find 
no evidence of correlations between cluster dynamical state and either the gas density slope in the inner regions or temperature, but 
do find some evidence of a correlation between dynamical state and outer gas density slope. We also find a weak correlation between 
dynamical state and both central gas normalisation and inner cooling times, but this is only significant at the 10% level. We conclude 
that, for the X-ray cluster population as a whole, both the central gas properties and the angle-averaged, large-scale gas properties are 
linked to the cluster dynamical state. We also investigate the central cooling times of the clusters. While the cooling times span a wide 
range, we find no evidence of a significant bimodality in the distributions of central density, density gradient, or cooling time. Finally, 
we present the gas mass-temperature relation for the REXCESS sample, finding that h{z)Mgas y'"*" which is consistent with 
the expectation of self-similar scaling modified by the presence of an entropy excess in the inner regions of the cluster and consistent 
with earlier work on relaxed cluster samples. We measure a logarithmic intrinsic scatter in this relation of ~ 10%, which should be a 
good measure of the intrinsic scatter in the Mg^^-T relation for the cluster population as a whole. 
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1. Introduction 

The X-ray emitting gas in galaxy clusters contains the signatures 
of important evolutionary processes such as mergers, AGN ac- 
tivity, galaxy interactions and tidal stripping. The impact of gas 
physical processes on the observable X-ray properties of galaxy 
clusters must be fully understood in order to use galaxy clus- 
ters to test the predictions of structure formation models, and to 
understand the relationship between AGN activity, galaxy evo- 
lution and the evolution of large-scale structure. The global scal- 
ing relations between cluster observables such as X-ray luminos- 
ity/temperature and cluster mass must be well constrained in the 
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local Universe so that cluster evolution to high redshifts can be 
investigated; this is not possible without a good understanding 
of the processes that lead to the observed scatter in these rela- 
tionships. 

Since the first evidence that the gas and dark matter content 
of galaxy clusters does not scale according to the simplest self- 
similar predictions (e.g. Edge & Stewart 1991; Arnaud & Evrard 
1999), non-gravitational heating processes have been thought to 
play an important role in determining the X-ray properties of 
galaxy clusters. Recent results have shown that the dark matter 
component of galaxy clusters does scale self-similarly, in good 
agreement with theoretical predictions (e.g. Pointecouteau et al. 
2005; Vikhlinin et al. 2006). It has also been shown that the de- 
parture of the observed Lx/Tx relation from theoretical predic- 
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tions is due to increased entropy in cluster gas (e.g. Ponman et al. 
2003; Pratt, Arnaud & Pointecouteau 2006). A combination of 
radiative cooling and feedback processes associated with galaxy 
formation, either from galaxy winds or AGN activity associ- 
ated with central supermassive black holes, is the most widely 
accepted explanation for these results (e.g. Voit 2005). Cluster 
mergers are also expected to have an important effect on their 
observed X-ray properties: for example, simulations have shown 
that cluster merger processes are likely to have an important ef- 
fect on the Lx/Tx relation (Rowley et al. 2004), although cooling 
and heating processes in the inner regions may be the dominant 
source of scatter (O'Hara et al. 2006). 

Radial surface brightness profiles and gas density profiles 
have been the key means of obtaining information about the 
structure and scaling properties of the intracluster medium 
(ICM) since the advent of X-ray imaging of clusters, providing 
an essential tool to enable study of three-dimensional gas dis- 
tributions, and, together with temperature, the study of the gas 
entropy distributions and the total mass profiles via the assump- 
tion of hydrostatic equilibrium. ROSAT studies (e.g. Neumann 
& Arnaud 1999) suggested that a /? model was an adequate 
parametrization of the X-ray surface brightness distribution out- 
side the core region in galaxy clusters, with clear evidence for 
self-similarity, but a large dispersion in the central regions. The 
core regions of clusters are now much better resolved with 
XMM-Newton and Chandra, and more complex models are re- 
quired to fit their surface brightness profiles (e.g. Pratt & Arnaud 
2002; Pointecouteau et al. 2004; VikhUnin et al. 2006). There is 
also some evidence that the slope of the gas density profile may 
steepen at large radii in some clusters (e.g. Viklinin et al. 1999; 
Neumann 2005; Vikhlinin et al. 2006). Constraints on the ra- 
dial distribution of hot gas in the galaxy cluster population are 
crucial both for accurate estimation of total cluster masses and 
for studies of gas entropy. Spherically symmetric methods for 
constraining the gas distributions of galaxy clusters have obvi- 
ous limitations when applied to samples that are not selected for 
regularity; however, as a straightforward method that is easy to 
apply to both observations and simulations, they remain an es- 
sential tool for studying cluster properties in three dimensions. 

In this paper we examine the gas distributions in 3 1 clusters 
drawn from the Representative XMM-Newton Cluster Structure 
Survey (REXCESS). Full details of the sample selection func- 
tion can be found in Bohringer et al. (2007). To achieve a statis- 
tically well-defined sample which fully samples the X-ray lumi- 
nosity function and is unbiased with respect to dynamical state, 
a sample of 33 clusters was constructed from the REFLEX cat- 
alogue (Bohringer et al. 2004) based on the following criteria: 
redshift z < 0.2; close to homogeneous coverage of luminos- 
ity space; kT > 2.0 keV; detectable with XMM-Newton to a ra- 
dius of ~ /?50(); and distances selected to optimise the extent of 
the cluster within the XMM-Newton field-of-view. In addition, 
a firm detection of more than 30 photons in the original RASS 
detection, and a low column density towards the source were 
also required. Observations of the full sample of 3 1 clusters plus 
2 archive observations have now been completed. The analysis 
presented here consists of the 31 clusters for which it is rea- 
sonable to carry out a 1 -dimensional analysis. We exclude the 
clusters RXCJ2152-1942 (Abell 2384B) and RXCJ0956-1004 
(Abell 901/902) which are highly uregular/multiple cluster sys- 
tems and hence cannot be analysed in this way. 

Throughout this paper we adopt a ACDM cosmology with 
Hq = 70 km S-' Mpc-', Qm = 0.3 and Qa = 0.7. 



2. Data analysis 

2.1. Data preparation 

Table [T] lists the global properties of the cluster sample and the 
details of the XMM-Newton observations. Observations were re- 
trieved and reprocessed with the XMM-Newton Science Analysis 
System (SAS) version 7.0, ensuring up to date calibration 
across the sample. Data sets were cleaned for flares, PATTERN- 
selected and corrected for vignetting as described in Pratt et al. 
(2007). The background used for the present analysis consists 
of custom event files generated from data taken in Filter Wheel 
Closed (FWC) mode, which correspond to an accurate repre- 
sentation of the particle and instrumental background present 
in XMM-Newton observations. These background data sets were 
cleaned, PATTERN-selected and vignetting corrected as above, 
then recast to have the same aspect as the observation data files. 
The background event lists were rescaled to the source quiescent 
count rate in the [10-12] and [12-14] keV range for EMOS and 
EPN cameras, respectively, by adjusting the WEIGHT column 
in each background event file. 

Sources other than the target object were detected in a broad 
band ([0.3-10.0] keV) coadded EPIC image using the SAS 
wavelet detection task ewavdetect, with a detection threshold 
set at 5o-. After visual screening, detected sources were excluded 
from the event file for all subsequent analysis. 

2.2. Surface brightness and gas density profiies 

Vignetting-corrected surface brightness profiles were extracted 
for each camera in the 0.3 - 2.0 keV band from source and 
scaled background event files in 3'.'3 bins out to a radius of 15'. 
The profiles were centred on the peak in X-ray surface bright- 
ness. The background surface brightness profiles were then sub- 
tracted, and the background subtracted profiles from the three 
cameras co-added. No weighting for camera response was ap- 
plied at this stage. We then applied the second stage of back- 
ground subtraction to account for the X-ray background. As the 
REXCESS sample was chosen with specific field-of-view crite- 
ria, all of the observations include an outer region which can be 
used to measure accurately the cosmic X-ray background (CXB) 
component. The residual background level due to the CXB in the 
co-added, background subtracted surface brightness profile was 
determined by a careful visual analysis of the region where the 
profile flattens due to background domination; the mean level in 
this region was then calculated and subtracted, and the resulting 
profile rebinned to a significance of 3cr per bin. 

Deprojected, PSF-corrected emission measure profiles were 
obtained from the surface brightness profiles using the non- 
parametric method described in Croston et al. (2006). The re- 
sponse matrices used were obtained using the Ghizzardi et al. 
(2001) parametrisation of the XMM-Newton PSF, weighted by 
the contribution of each camera to the combined profile. The pro- 
files were then converted to gas density by calculating a global 
conversion factor for each profile in XSPEC using the global 
temperatures listed in Table[T] which are spectroscopic tempera- 
tures estimated in the [0.15 < < 1] Rsm aperture. A correction 
factor to take into account radial variations of temperature and 
abundance was calculated for each radial bin using a parameter- 
isation of the projected temperature and abundance profiles, as 
detailed in Pratt & Arnaud (200313. The surface brightness pro- 
files and corresponding gas density profiles for each cluster are 

' The use of projected values here is acceptable, as the temperature 
dependence of emissivity in the 0.3 - 2.0 keV band is < 5%, and there 
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Table 1. Cluster properties. Columns: (1) Cluster ID, (2) redshift, (3) column density in 10^° cm"^, (4) X-ray temperature 
from a mekal fit in the [0.15 - 0.75]/?5oo region in keV, (5) X-ray temperature from a mekal fit in the [0.15 - l]^5oo region 
in keV, (6) abundance obtained from a mekal fit in the [0.15 - l]/?5oo region, in units of solar abundance, (7) Rsqq in kpc, 
iteratively determined from the M500 - Yx relation of Arnaud, Pointecouteau & Pratt (2007) as described in the text, (8) 
logarithm of gas mass to TJgoo determined as described in the text in Mq, (9) inner slope of the gas density profile determined 
in the region with r < 0.057?5oo, (10) outer slope of the gas density profile determined in the region[0.3 - 0.8]7?5oo as defined 
in the text. 
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presented in the Appendix. The gas density profiles in tabular 
form are included in the electronic version of this paper. 

3. Results 

3.1. Global properties of cluster gas density profiles 

Fig. [U shows the 31 cluster gas density profiles superposed to 
allow an investigation of their global properties. The individual 
surface brightness profiles (in the energy range 0.3 - 2.0 keV) 
and gas density profiles are given in the Appendix. There is a 
considerable amount of scatter in the unsealed profiles (top left 
panel), which is unsurprising. Structure formation models pre- 
dict that the mean cluster density scales with pdz) and so ac- 
cording to h{z)^ and predict a universal profile shape when the 
radial coordinate is scaled according to the cluster mass. We use 
a scaling radius of R5Q0, defined as the radius enclosing a mean 
overdensity of 500 times the critical density. We therefore expect 
that gas density profiles scaled by h(zy^, and in radial units of 
r/7?5oo should coincide. 

The physical distances were converted to scaled radius us- 
ing the values of Rsoo given in Table [T] which were determined 

are no significant differences between projected and deprojected abun- 
dance. 



iteratively as described by Kravtsov et al. (2006), based on the 
A/500 - Yx relation of Arnaud, Pointecouteau & Pratt (2007) (see 
also Maughan et al. 2007). Here Yx is the product of the gas mass 
within 7?5oo and the temperature in the [0.15- l]7?5oo region. This 
approach allows us to determine R^m directly from the data with- 
out first fitting total mass profiles. For self-consistency, we prefer 
to use a scaling relation determined from XMM-Newton results, 
rather than relying on results determined from lensing masses, 
which have large individual uncertainties and are only applica- 
ble to the high mass regime, or from simulations that do not yet 
reproduce all of the observed properties of the cluster popula- 
tion. A comparison between the empirical relation of Arnaud et 
al. (2007) and the simulated relation of Nagai et al. (2007) shows 
only a small normalisation offset (< 8%), so that the only differ- 
ence to our results if the simulated relation were used would be 
to shift all the profiles in the direction of smaller radii. It would 
therefore have no effect on any of our conclusions, which are all 
based on the relative profiles. 

The total gas mass to a radius of Rsm was determined for the 
sample by integrating over the gas density profiles. For profiles 
where the data does not extend to /?5oo (9/3 1 clusters) we extrap- 
olated the profiles from the outer bin of the density profile (never 
less than 0.8/?5oo) assuming a power law slope in log space, de- 
termined by a fit to the outer four bins. The errors on gas mass 
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were determined by combining in quadrature the errors obtained 
from Monte Carlo simulation over the statistical errors in the 
surface brightness profiles and the errors due to the temperature 
uncertainties. 

The scatter in gas density at 0.3/^500 is crlirie) - 0.265. We 
then scaled the gas densities according to their expected evolu- 
tion with redshift: oc m)h{z)^ (shown in the middle left panel of 
Fig. [Til, which reduced the scatter at 0.3/?5oo to cr/{ne} = 0.236. 
The bottom right panel of Fig. [T] shows the relative dispersion as 
a function of radius for this scaling (in red), which drops from 
~ 100% in the inner regions to a minimum of ~ 15% at 0.6/?50(), 
before increasing slightly in the outer regions. 
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3.2. Dependence of gas density on temperature 

Fig.|2]shows the evolution-corrected gas density profiles colour- 
coded by cluster temperature. At the smallest radii probed by 
this data, the dispersion in profiles is uncorrelated with temper- 
ature (e.g. left hand panel of Fig. O. We investigate the origin 
of this scatter in §4.3. At intermediate radii (~ 0.3R^qo), there 
is a clear systematic bias with system temperature, although this 
trend becomes weaker again at the maximum radii probed by the 
data. Fig. [3] also quantifies the relationship between gas density 
normalisation at a radius of 0.3/?5oo and at 0.7/?5oo. Results are 
shown for the redshift-scaled gas density profiles; however, us- 
ing the unsealed profiles does not significantly alter the result. 
There is evidence for a correlation in both cases, although the 
trend is weaker at Q.lRsoo- We obtained a null hypothesis prob- 
ability < 0.0001 on a Spearman rank test for the comparison 
at 0.3 /?50() and a higher null hypothesis probability of 0.004 at 
0.7/?5()o. We carried out orthogonal linear regression, finding a 
slope of 0.50 + 0.08 for the relationship between ne(Q.3 /?5oo) 
and T. 

The departure of the density from the simple scaling is best 
discussed in terms of the radial entropy profiles of the sys- 
tems. To be consistent with previous work, we define entropy 
as 5 = T/n'^J . The observed departure from simple scaling of 
the entropy-temperature relation for clusters (e.g. Finoguenov et 
al. 2002, Ponman et al. 2003) implies a temperature dependence 
of the gas density normalisation such that oc t'O-SZs jj^jg jg jjj 
excellent agreement with the trend we see at 0.3 Rsoq. However, 
at 0.7 R5Q0, the correlation is weaker and we find a flatter slope of 
0.25 + 0.06. It is apparent from the righthand panel of Fig. [3]that 
the correlation at 0.77?5(K) is mainly due to lower densities in the 
systems below kT ~ 3.0 keV; at higher temperatures the relation 
appears to be flat. This suggests that the "entropy excess" may 
extend to larger scaled radii in cooler systems, consistent with 
the expectation that non-gravitational processes have a greater 
efl'ect at the low mass end of the cluster population. 

We then scaled the density profiles by T'^-^^^ (in addition to 
the redshift scaling discussed above), as shown in the middle 
right panel of Fig. [1] As shown in the bottom right panel, the rel- 
ative dispersion decreases for this scaling (shown in red) out to 
a radius of ~ 0.57?5()o, reaching a minimum of cr/{ne} = 0.133 at 
0.57?5oo. Beyond ~ O.SRsm the relative dispersion for this scaling 
is slightly higher than for the redshift-evolution-only profiles. 
This result is consistent with the weaker correlation between gas 
density and temperature at 0.7/?5oo, supporting the conclusion 
that an entropy excess is likely to be significant in the cluster 
inner regions only. 
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Fig. 2. Gas density distributed scaled according to expected evo- 
lution with redshift, and coded by temperature so that blue cor- 
responds to a temperature of ~ 2.0 keV and red a temperature of 
~ 8 keV. 

3.3. Radial dependence of density profile slope 

We investigated the radial dependence of the density profile 
slope, defined as the logarithmic gradient of density with ra- 
dius, and determined from the deprojected profiles, with a 
binning chosen to correspond to a typical temperature profile 
for the sample. In Fig. HI we present the superposed plots of 
dlog(ne)/dlog(r) for the cluster sample, colour-coded by temper- 
ature in the same way as for Fig.|2] There is a large scatter at all 
radii of the profiles, which in many cases are not very smooth. 
The scatter at a radius of 0.3/?5oo is cr/ (a) - 0.20. It is evident 
from Fig. |4] that in most cases a appears still to be increasing in 
the outermost regions (e.g. at a distance of ~ 0.8 - l.O/^soo) in 
contrast to the power-law behaviour expected for a y6-model pro- 
file. This is in agreement with earlier work (e.g. Vikhlinin et al. 
1999; Neumann 2005), and has important implications for clus- 
ter mass estimates obtained using analytical models. We have 
fully taken into account systematic uncertainties in the back- 
ground level in determining the outer slopes: the systematic un- 
certainty in the background level is typically ~ 2 - 3%, leading 
to an uncertainty in outer slope that is negligible compared to 
the statistical uncertainties on surface brightness. In one case, 
RXCJ2234.5-3744, the background uncertainty is ~ 14% due 
to residual flare contamination, which leads to a ~ 3% error in 
the [0.3 - 0.8]7?5oo slope. In FigH] the mean slope is shown only 
out to 0.97?g(K), where the effects of systematic uncertainty in the 
background level are negligible. Our conclusions relating to the 
gas density gradient are therefore robust. 

3.4. Parametrization of cluster structure 

In order to further investigate the distribution of cluster struc- 
ture and the dependence of cluster structure on other properties 
such as gas temperature and dynamical state, we first consid- 
ered the outer logarithmic slope of the gas density profiles. We 
chose to use the region of [0.3 - 0.8]/?5oo, since, as shown in 
Fig. [1] and Fig. |4] all of the profiles have steepened from their 
inner slopes by 0.37?5oo, and all extend to at least 0.87?g()o. The 
use of a fixed outer radius minimizes the possibility of bias due 
to brighter, and hence higher temperature, clusters extending to 
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Fig. 1. Gas density profiles for the entire sample. Top left: unsealed profiles; top right: profiles of unsealed density with scaled 
radius; middle left: profiles of density scaled for redshift evolution; middle right: profiles of density scaled according to r^^^s 
implied by modified self-similar S-T scaling; bottom left: density profiles with a representation of the 1-sigma scatter of the sample; 
bottom right: the relative dispersion (cr/ {11^}) as a function of radius for the cluster sample with profiles scaled according to expected 
evolution with redshift (black) and by T'^-^^^ (red). 
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Fig. 3. Plots of gas density normalisation vs. cluster temperature at 0.03/?5oo (left), 0.3/?5oo (middle) and at O.lRsoo (right), for the 
redshift-scaled density profiles. For the latter two cases, lines of best fit obtained as described in the text are overplotted. 
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Fig. 4. Profiles of a, the density profile slope for the entire cluster 
sample, colour-coded by cluster temperature as in Fig.|2] Thick 
red line indicates the mean profile. 
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Fig. 6. Comparison of the inner and outer slopes of the density 
profile (a<o.o5 and ao3-o.&, respectively) for each cluster 



larger radii. We express the slope as a y6 value, where /3 is defined 
as 1/3 times the slope, for comparison with other works using 
analytical models. The outer slope in the region [0.3 - 0.8]/?5()o 



is hereafter referred to as ao.3-o.8> or y6o.3-o.8 when expressed as 
a y6 value. We also measured the inner logarithmic slope within 
0.05/?5oo, before any of the profiles have begun to steepen signif- 
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icantly, which is hereafter referred to as Q'<().()5. Both ao.s-o.s and 
Q'<o.05 were measured by a simple linear least-squares fit in log 
space over the relevant region. As discussed above, the effects of 
systematic uncertainties in the background level have been fully 
taken into account and are less than the statistical uncertainties 
in nearly all cases. 

In Fig. |5] (top) we plot the distribution of j0o.3-o.8 and q'<o.o5 
for the gas density profiles scaled according to expected evolu- 
tion with redshift. The mean value for q'<o.o5 is -0.5 ± 0.3, with 
the large scatter reflecting the range from clusters with central 
cores to those that are highly peaked (using a fixed inner bound- 
ary in scaled radius does not reduce the scatter or significantly 
alter the mean slope). We are confident that the values of aoj-o.s 
are unaffected by PSF-correction or resolution issues: in Croston 
et al. (2006) we showed that the PSF correction is very accu- 
rate beyond ~ 0.05 arcmin, and we find that Q'0.3-0.8 is uncorre- 
cted with redshift or angular scale. The mean value of aoj-o.s 
is -1.80 ± 0.28 (corresponding to y8().3_o.8 - 0.60 + 0.10), which 
is consistent with other studies, e.g. Neumann & Arnaud (1999), 
Ota & Mitsuda (2004), Vikhlinin et al. (2006). The scatter in 
Q'0.3-0.8 is considerably smaller than for Q'<o.o5; however, a few 
clusters are exceptionally flat in their outer regions (three clus- 
ters have /3 < 0.5). There is no evidence for bimodality in a<o.o5- 
We also investigated whether q'<o.o5 and Q'0.3-0.8, the inner and 
outer slopes, are coiTelated. As shown in Fig. |6l no such corre- 
lation appears to exist, indicating that in general the dynamical 
evolution of the angle-averaged gas distribution in the central 
and outer regions of the clusters are not closely connected. 

The REXCESS sample includes clusters with evidence for 
irregularity, and in some cases the choice of the centre for the 
profiles is very dependent on the radius used for centroiding 
(see Section |4~2l ). We therefore examined the effect of this uncer- 
tainty on the structural parameters of the density profiles by con- 
sidering three of the most irregular clusters, RXCJ2048-1750, 
RXCJ2129-5048 and RXCJ1516-0056. For each of these clus- 
ters, surface brightness profiles were extracted using centroid po- 
sitions obtained within radii of 0.3, 0.5 and 1 .0 times R5Q0, which 
were deprojected in the same way as for the original profiles. We 
found that the choice of centroid unsurprisingly has an important 
effect on the measured value of Q'<o.o5, introducing a scatter of 
between 40 - 67 percent; however, the results for the inner two 
choices of centroiding region, more appropriate for studying the 
slope in the inner regions, were in reasonable agreement. The 
choice of centroid does not appear to introduce significant uncer- 
tainty in Q'0.3-0.8: the dispersion between the three cases ranged 
from 2-6 percent in these "worst-case" clusters. We comment 
below on the effect of centroid choice on other results. 



3.5. Dependence of cluster structure on temperature 

In Fig. |2]we plot the inner and outer density slopes (q'<o.o5 and 
P03-0.&) against the cluster temperature in the region [0.15 - 
l]7?50o- /So.3-0.8 appears to coiTelate with global temperature. 
Based on a Spearman rank test, we find a ~ 2 percent proba- 
bility that such a coiTelation could occur by chance. There is no 
significant correlation between a<o.o5 and the global cluster tem- 
perature. It is clear from Fig. |4] that the slope does not depend 
on temperature over most of the radial range. It is only at radii 
beyond ~ 0.3/?5()o that the slopes of the cooler profiles are pref- 
erentially below the mean, with hotter profiles mainly above the 
mean. 
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Fig. 8. Cluster gas mass vs. global temperature. Solid line is the 
best fitting relation to this sample obtained using WLSS, dashed 
line is the results obtained with BCES, and the dotted line is the 
best-fitting relation for the sample of regular clusters discussed 
in Arnaud et al. 2007. 



4. Discussion 

4. 7 . Total gas mass and the Mgas-T relation 

Fig.[8]shows the relationship between gas mass and global tem- 
perature for the REXCESS sample. Here we use temperature 
in the [0.15 - 0.75]/?5oo region to enable direct comparisons 
with previous work. We fitted a scaling relation of the form 
h{z)Mgas - C[Txl5keVy using two linear regression methods, 
a weighted least squares method that incorporates an intrinsic 
scatter term in the errors (WLSS), and BCES, which is an or- 
thogonal method that does not take into account the statistical 
errors. These methods are discussed in more detail in Pratt et al. 
(2006). The results of the fits are given in Table |2l The slopes 
and normalisations obtained with the two methods are consis- 
tent, and are also in good agreement with those obtained for a 
smaller sample of relaxed clusters by Arnaud et al. (2007). The 
results of our fits and the Arnaud et al. relaxed cluster relation 
are shown in Fig[8] As shown in Table|2] the intrinsic logarithmic 
scatter for the REXCESS sample is a factor of 2.5 times higher 
than for the Arnaud et al. (2007) relaxed cluster sample. This 
is likely to arise as a result of the wider range of cluster mor- 
phologies in our sample, although the sample appears to include 
clusters with both higher and lower gas mass for their tempera- 
ture. 

We investigated the effect of choice of cluster centre on the 
gas mass measurements and found that for the three most dis- 
turbed systems the gas masses obtained varied by between 7 
and 12% for centroiding using regions of radius 0.1, 0.3 and 0.5 
x/?5oo. As the cluster dynamical state appears to be relatively 
independent of temperature (see Section |4!2j, we do not expect 
this systematic uncertainty to aff'ect our conclusions about the 
scaling of gas mass and temperature. As, in addition, our sample 
selection strategy is unbiased with respect to cluster dynamical 
state, the ~ 9% scatter about the Mgas-T relation we measure for 
this sample should be a good measure of the intrinsic scatter in 
the relation for the cluster population as a whole. 
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Fig. 7. Plots of the inner (left) and outer (right) slopes (q;<o.o5 and /Joj-o.s), respectively, versus the global cluster temperature, Tx- 

Table 2. Results of linear regression analysis for the Mgas-T relation, y is the slope in log-log space, C the normaUsation, and cTrat 
and (Tuitrms the raw and intrinsic logarithmic scatter about the relation, respectively. 



Sample 


method 


r 


log{C) 






REXCESS 


WLSS 


1.986±0.111 


13.652 ±0.020 


0.0928 


0.0903 




BCES 


2.122 ±0.121 


13.661 ±0.019 


0.0989 


0.0962 


Arnaud et al. (2007) 


WLS 


2.10 ±0.05 


13.65 ±0.01 


0.048 


0.036 



The steep relation between gas mass and temperature im- 
pUes a decrease in the total gas mass content of cooler clusters 
relative to higher mass systems, when compared with standard 
predictions. This is clearly connected with the excess of entropy 
seen at intermediate radii in cooler systems and their relatively 
flat outer slopes. All of these factors suggest that significant non- 
gravitational heating is likely to have occurred, raising the gas 
entropy and lifting material beyond Rsm. Energetically it is most 
appealing if this energy injection occurred prior to cluster col- 
lapse (eg., McCarthy et al 2007), but a single level of "preheat- 
ing" is unable to explain the full range of observed properties, 
both as a function of cluster temperature and as a function of 
the variations within in a temperture band. We will explore the 
theoretical models that can account for the observed trends in a 
future paper. 



sity slope (a<o.o5); however, there is evidence for a correlation 
between Pi/Po and y6o.3-().8^ with a null hypothesis probability of 
~ 0.1% on the generalized Kendall's r test. There is weak evi- 
dence for a correlation between P3/P0 and temperature, with a 
null hypothesis probability of ~ 10% on a generalized Kendall's 
T test. This correlation likely arises as a result of the correlation 
between /? and temperature, as discussion in Section 13.51 The 
relations between gas structure and P3/P0 are shown in Fig.|9] 



4.2. Dependence of cluster structure on dynamical state 

The radial gas distribution in galaxy clusters is likely to be 
strongly aff'ected by the cluster dynamical state, since mergers 
are expected to disrupt the gas structure significantly. We used 
two different quantitative measures of substructure in the 2-D 
surface brightness distribution as means of characterising the 
cluster dynamical state: the power ratio method of Buote & Tsai 
(1995) and centre shifts (e.g. Mohr et al. 1995, Poole et al. 2006). 

The power ratio method used here is described in more de- 
tail in Pratt et al. (2007). Here, we make use of power ratio mea- 
surements for the entire REXCESS sample (to be discussed in 
a forthcoming paper), determined in an aperture of radius /?5oo- 
We examined the dependence of cluster structure on the three 
lowest order power ratios of relevance for this work: P2/P0, 
which corresponds to a measure of ellipticity, P^/Po, which is 
the best measure of further substructure, and P4/P0, which also 
measures deviation from a relaxed dynamical state. As some of 
the Pi/P[) values are formally upper limits, we used the gener- 
alized Kendall's t test for censored data, as implemented in the 
ASURV package. We did not find a correlation with inner gas den- 



Centre shifts were determined for the entire REXCESS sam- 
ple. Centroids were obtained for regions of « x 0.1 x R^m with 
n - 1..10 (i.e. the region with r < O.lRsm was not included, 
so that the centre shifts are not affected by the large scatter in 
central gas properties). The results of this analysis will be pre- 
sented in more detail in a forthcoming paper; however, for the 
purposes of a comparison with gas structural properties we used 
(w), defined as the standard deviation of the projected separa- 
tions between the X-ray peak and centroid at each radius (e.g. 
Maughan et al. 2007) in the region between 0.1 x R^oq and TJgoo. 
We found that the exclusion of the central region, possibly af- 
fected by a cooling core, did not significantly affect the measured 
(w) values, except in the cases where (w) is very small. In Fig.|9] 
we plot (w) against inner and outer gas density slopes (q;<o.o5 
and yS().3_o.8)- On a Spearman rank test, we find no strong evi- 
dence for a significant correlation with Q'<o.o5, but find evidence 
for a correlation between (w) and j0o.3-o.8, with a null hypoth- 
esis probabilty of 0.08% on a Spearman rank test, (w) is also 
uncorrected with temperature. The lack of correlation with in- 
ner density slope suggests that cluster substructure and central 
cooling behaviour are independent. 



J.H. Croston et al.: REXCESS gas density distributions 



01 

in 
o 



V 



o 



-9 -8 -7 
log(P3/Po) 



d 



^ CO 

O 
in 

ro 

d 

A "3 

in 
d 



d 



<i- <i- 



-8 -7 
iog(P3/Po) 



o 



> 

0) CD 



-8 -7 
log(P3/Po) 



o 
d 



+ 



-2 -1 
log(<w>) 



03 

d 
d 



o 

A «3 



d 



.•■t-.% •) 

+ J* 

J , , L_ 

-2 -1 
log(<w>) 



> oo 

(D 



o 



♦ + 

•t: • 

. • t 

■ * 

J 1 , 1 , I 

-2 -1 

log(<w>) 



Fig. 9. 



The relationship between cluster dynamical state as parametrised by power ratios and centre shifts and gas properties. Top: 
bottom: (w). Left to right: inner slope of gas density (q'<o.o5), outer slope of gas density (jSoj-o.s), temperature. 
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Fig. 10. Cooling time profiles for the cluster sample. Left: unsealed profiles in physical units; middle: unsealed gas density in units 
of scaled radius (corresponding to the top left panel of Fig. [1]); right: profiles obtained from the redshift-scaled gas density profiles 
(corresponding to the middle left panel of Fig. [TJ. 



4.3. Cooling times 

Cooling time profiles were determined from the gas density and 
temperature profiles for each cluster, where the cooling time, 
tcooi is defined as: 

'-"'^'^ = — !;;(?) — 

where Lx{r) is the X-ray luminosity at radius r determined us- 
ing the appropriate mekal model parameters at each radius in 
xsPEC. Fig. [To] shows the cooling time profiles for all 31 clus- 



ters in physical and scaled units in the radial range where the 
temperature profile is well constrained for each cluster. As the 
temperature profiles we use here are not deprojected, we com- 
pared the cooling time profiles obtained using projected and de- 
projected temperature profiles for the cluster with the steepest 
central temperature drop. We conclude that the use of projected 
temperatures introduces an uncertainty of at most ~ 15% at radii 
less then Q.QlRsm and has negligible effects at larger radii. 

The lower envelope of the cooling time profile distribution 
for REXCESS is in good agreement with previous work (e.g. 
Sanderson et al. 2006, Bauer et al. 2005); however, the dis- 
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Fig. 11. Histograms of central gas density and central cooling time (at 0.03/?5()o, the innermost radius at which both density and 
temperature data is available for all clusters), showing that there is no strong evidence for bimodality in their cooling properties. 
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Fig. 12. Comparisons between central gas properties and dynamical state. Left: P3/P0 vs. ne(0.007/J5oo) (top) and frao/(0.03/?5oo) 
(bottom); right: <w> vs. riei^miR^^)) (top) and fcoo/(0.037;goo) (bottom). 



persion in the inner regions of this sample is higher than seen 
in other samples, particularly those dominated by cooling core 
clusters (e.g. Voigt & Fabian 2004, Sanderson et al. 2006), with 
a value of crl{tcooi) - 0.85 at 0.037?5oo (as the cooling time pro- 
files are dominated by the behaviour of the gas density profiles, 
the scatter as a function of radius is similar to that shown in 



Fig. [TJ. We find a mean slope in the region [0.2 - 0.8]7?5oo of 
1.4 + 0.3, which is consistent with the results of Bauer et al. 
(2005). In the region [0.05 - 0.2]/?5oo the mean slope for our 
sample is ~ 0.9 + 0.3, somewhat flatter than that seen over a sim- 
ilar range in scaled radius by Sanderson et al. (2006), due to the 
larger scatter towards higher cooling times in our sample. The 
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Fig. 14. The gas density normalisation (at 0.37?5oo) vs. temperature distribution for the simulated cluster sample using both the 
spectroscopic-like temperature (left) and the emission-weighted temperature (right). No significant correlation is seen 
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Fig. 15. The outer slope (jSoj-o.s) vs. temperature distribution for the observed and simulated cluster samples using both the 
spectroscopic-like temperature (left) and the emission-weighted temperature (right) for the simulated sample. 



spatial resolution of our temperature profiles in the inner regions 
limits our ability to draw strong conclusions about the central 
cooling properties of this sample; however, the large dispersion 
in central cooling properties suggests that there are more clus- 
ters with high central cooling times in this sample compared to 
other studies. Nevertheless, tcooi drops below the Hubble time at 
a radius of ~ 50 - 200 kpc for virtually all of the clusters in the 
sample. 

Traditionally, clusters have been divided into two classes ac- 
cording to their central cooling times: "cool-core" systems with 
short central cooling times and high central densities and "non- 
cool-core systems" with cooling times comparable to the Hubble 
time. Because of their higher surface brightness, cool-core sys- 
tems have tended to be over represented in previous deep ob- 
servations of clusters, leading to the perception that the distribu- 
tion of system cooling times is bimodal. However, as shown in 
Fig. [TT] the distributions of gas densities and cooling times for 
the REXCESS sample at the innermost radius appear to be well 
represented by a broad single peaked distribution, and any divi- 
sion of this sample based on the central cooling time would be 
arbitrary. We note, however, that our cooling time profiles do not 
extend to very small radii, and so we cannot rule out bimodal- 
ity in the cooling time profiles on smaller scales. In particular 
there does appear to be a subset of profiles which continue to 
have a steep gradient in cooling time at small radii, which corre- 



sponds primarily to those profiles with fcoo/(0.03/?5oo) < 1 Gyr. 
The mean slope for profiles with inner cooling times higher than 
1 Gyr is significantly flatter, with several profiles that are asymp- 
totically flat in the central regions. Hence our data do not rule out 
models in which conduction can stabilize the cores of clusters 
with high central cooling times (e.g. Donahue et al. 2005). As 
noted in Section [331 the distribution of inner gas density slope 
(Q'<o.05) also shows no strong evidence for bimodality. These 
conclusions are independent of the choice of radius and gas den- 
sity scaling. 

Finally, we investigated the connection between central gas 
properties and cluster dynamical state by comparing the inner 
cooling times and gas densities with (w) and P^IPq. Although it 
is widely anticipated that clusters with long central cooling times 
result from cluster mergers, recent papers have suggested that 
central density cusps are unlikely to be destroyed in the merger 
event (eg., McCarthy et al. 2004, Borgani et al. 2004, Poole et 
al., 2006). It is therefore of considerable interest to investigate 
the connection between central gas properties and cluster dy- 
namical state by comparing the inner cooling times and gas den- 
sities with (w) and P^/Po- Fig.[T2]illustrates that coiTelations are 
present, in the sense that clusters with lower central densities, 
and higher central cooling times, tend to show more evidence 
for disturbance. However, the trend is driven by a few systems 
with the highest level of disturbance, and the null hypothesis that 
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Fig. 13. Comparison of the observed gas density profiles ( 1 cr dis- 
persion in red) and simulated profiles of Borgani et al. (2004) 
(dashed lines) scaled by p,,,>- 



there is no coiTelation can only be rejected at the ~ 85 - 97% 
level using the Spearman rank correlation test (for {w}) or the 
generalized Kendall's r test (for Pj,/Pq). While there is there- 
fore weak evidence that merger activity can affect cluster central 
densities and cooling times, it does not appear to affect the slope 
of the density profile (Section [4.21 ) in the inner regions. Clearly 
these issues need to be addressed by careful comparison with 
cosmological simulations. 

4.4. Comparison with simulations 

We compared the observed gas density profiles with the simu- 
lated profiles of all clusters with kT > 2 keV from Borgani et 
al. (2004), in which the SPH code GADGET-2 (Springel 2005) 
was used to simulate a concordance ACDM model (Q^ - 0.3, 
Qa = 0.7, o-g = 0.8, h = 0.7) within a box of 192 h'^ Mpc 
on a side, using 480^ dark matter particles and an equal num- 
ber of gas particles. The simulation included radiative cooling, 
star formation and galactic ejecta powered by supernova feed- 
back. Fig.[T3]shows in red the Icr dispersion of the entire sample 
of observed profiles as in the bottom left panel of Fig. [1] with 
the simulated profiles as dashed lines. Observed profiles were 
scaled by pi,,>, defined as lOOpcr/r.o- The agreement is very good 
in the radial range [0.02 - 0.3]/?5oo; however at larger radii the 
observed profiles are slightly flatter than the simulated profiles. 
We found that the simulated profiles have a mean slope in the 
radial range [0.3 - 0.8]7?5oo of -1.97 + 0.25, which is consistent 
within the errors with our measured mean slope of -1.80 + 0.28, 
but is slightly higher. For the purpose of this comparison we used 
the "true" R5Q0 values from the simulations. The small difference 
between the slope of the observed and simulated profiles at large 
radius may be due to differences in the definition of /?50(). We 
cannot apply our method for determining R^qo based on the ob- 
served Yx - M500 relation directly to the simulations, because 
spectroscopic-like temperature measurements that exclude the 
central regions of the cluster are not available. The relative dis- 
persion in the simulated profiles at 0.3/?5oo is (T/(«g) = 0.13, 
which is significantly lower than the value of ~ 0.22 for the 
redshift-scaled observed profiles. 



Nagai et al. (2007) carried out a similar comparison of the 
gas distributions in simulated and observed clusters, using an 
observed sample consisting of the z - relaxed cluster sample 
of Vikhlinin et al. (2006). They found slightly better agreement 
in the outer regions between the Chandra profiles and their sim- 
ulations that included cooling and star formation. Their slightly 
better agreement between observations and simulations may be 
a result of comparing with a sample including only relaxed clus- 
ters. Differences in the numerical treatments and implementa- 
tions of cooling and star formation in the simulations may also 
be relevant. 

We also investigated whether the temperature-dependence of 
gas density normalisation seen in the REXCESS clusters (see 
Section [32J is also present in the simulations. Figure [14] shows 
the relationship between gas density normalisation and temper- 
ature (shown for both the emission- weighted and spectroscopic- 
like temperatures) for the simulations. There is no significant 
temperature-dependence of gas density normalisation for the 
simulated cluster sample. 

Finally we compared the outer slope with gas tempera- 
ture, using both the emission-weighted temperatures within TJgoo 
and the spectroscopic-like temperatures defined in Rasia et al. 
(2005). Fig. [15] shows that the correlation between gas density 
slope O0O.3-O.8) and temperature observed in the REXCESS sam- 
ple is not present in the simulated data (on a Spearman rank test 
we find null hypothesis probabilities of 35% and 75% for T,/ and 
Ten; respectively, compared with < 1 % for the observed sample) 
- there are clearly a number of simulated clusters in the low T, 
steep y6 region of Fig. [15] which is not populated by the observed 
sample. However, the temperature measurements for the simu- 
lated sample do not exclude the central region which makes it 
difficult to draw firm conclusions. It is also clear from Fig. [15] 
that the temperature distribution in the simulated sample is dif- 
ferent to that of the observed sample, with fewer high tempera- 
ture clusters. 

These differences hint that the simulations may not match 
the true thermal history of the intracluster medium. In particu- 
lar, the heat input appears too centrally concentrated so that the 
excess entropy seen in the lower temperature clusters at interme- 
diate radii is not reproduced. A more detailed comparison with 
simulations will form part of a later paper. 

5. Conclusions 

We have presented the first detailed study of the structural prop- 
erties of cluster gas in a large, representative sample of nearby 
galaxy clusters. As the sample was selected by X-ray luminos- 
ity, it includes clusters of all dynamical states, allowing us to 
investigate the effect that the inclusion of less regular systems 
has on results obtained with previous studies of regular clusters. 
We found the following results: 

- The 1-D gas density profiles scale self-similarly, with a scat- 
ter ranging from ~ 100% in the inner regions to ~ 20% at 
a radius of 0.3/?5oo when expected evolution with redshift is 
taken into account. 

- Gas density normalisation at 0.3/?5oo is correlated with 
global cluster temperature, with a scaling of ne(0.37?g(K)) oc 
T'^^, consistent with the expectation of modified entropy- 
temperature scaling models. Using this scaling reduced the 
scatter in the gas density profiles at 0.3/?5()o to ~ 15%; how- 
ever, the scatter at larger radii is slightly increased, which 
indicates that the entropy excess is much less significant be- 
yond ~ 0.5/?5oo- 
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- The gas density slope continues to increase with radius in 
the region 0.5/?5()o - R500, as found by others, which is of 
importance for cluster mass estimates. 

- The outer gas density slope is correlated with X-ray temper- 
ature, primarily due to a lack of hot clusters with flat gas 
distributions. The flatter slope in lower temperature systems, 
combined with the entropy excess at intermediate radii and 
the steep slope of the Mgas-T relation suggests that gas has 
been displaced from the centres of the lower temperature sys- 
tems to larger radii. 

- Based on a characterisation of the cluster dynamical state us- 
ing power ratios and centre shifts, there is evidence of a cor- 
relation between cluster dynamical state and outer gas den- 
sity slope, and no correlation with inner gas density slope. 
There is also evidence for a correlation between dynamical 
state and the central gas properties (gas density normalisa- 
tion at 0.0077?g(K) and cooling time at 0.03/?5oo)- 

- There is no evidence for bimodality in the central gas den- 
sity, gas density slope or cooling times for this sample, sug- 
gesting that X-ray clusters form a single population with a 
continuous distribution of central gas properties. 

- The gas-mass temperature relation for the REXCESS sample 
is in good agreement with predictions of self-similar scaling 
modified by the presence of an entropy excess, and with pre- 
vious work on samples of regular clusters; however, the in- 
trinsic scatter is a factor of ~ 2.5 times higher than for the 
relaxed cluster population. 

- The scaling properties of the gas density profiles appear to be 
in broad agreement with those of simulated cluster samples 
(Borgani et al. 2004) at intermediate radius, with a slightly 
flatter slope in the outer regions. However, in contrast to the 
observational data there is no correlation between gas den- 
sity normalisation and temperature in the simulated sample, 
or between the outer slope and temperature. These discrep- 
ancies suggest that the non-gravitational heating of the intra- 
cluster medium may be too centrally concentrated in these 
models. 
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Fig. 16. Co-added MOSl, MOS2 and pn surface brightness profiles for the entire sample in the energy band 0.3 - 2.0 keV. 
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Fig. 17. Co-added MOSl, MOS2 and pn surface brightness profiles for the entire sample (cont.) 
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Fig. 18. Co-added MOSl, MOS2 and pn surface brightness profiles for the entire sample (cont.) 
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Fig. 19. Individual gas density profiles for the entire sample. Vertical lines indicate R500 for each cluster. 
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Fig. 20. Individual gas density profiles for the entire sample (cont.) 
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Fig. 21. Individual gas density profiles for the entire sample (cont.) 



